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Abstract. We numerically study the phase structure of two types of triangulated spherical surface models, 
which includes an in-plane shear energy in the Hamiltonian, and we found that the phase structure of the 
models is considerably influenced by the presence of the in-plane shear elasticity. The models undergo a 
first-order collapsing transition and a first-order (or second-order) transition of surface fluctuations; the 
latter transition was reported to be of second-order in the first model without the in-plane shear energy. 
This leads us to conclude that the in-plane elasticity strengthens the transition of surface fluctuations. 
We also found that the in-plane elasticity decreases the variety of phases in the second model without the 
in-plane energy. The Hamiltonian of the first model is given by a linear combination of the Gaussian bond 
potential, a one-dimensional bending energy, and the in-plane shear energy. The second model is obtained 
from the first model by replacing the Gaussian bond potential with the Nambu-Goto potential, which is 
defined by the summation over the area of triangles. 

PACS. 64.60.-i General studies of phase transitions - 68.60.-p Physical properties of thin films, nonelec- 
tronic - 87.16.D- Membranes, bilayers, and vesicles 



■ 1 Introduction 

I Surface models are conventionally defined by a surface 
. tension energy and a bending energy [Tl[51[3] . Both of the 
' energies play a role of maintaining the shape of surface 
against environmental external forces including thermal 
fluctuations (4ji5ji6ji7ji8,i£| , Thus, the surface models are 
always constructed to have resistance against tensile de- 
. formations and bending deformations, while no resistance 
' is assumed against in-plane shear deformations. This final 
assumption seems valid in fluid membranes at least, be- 
cause the in-plane shear deformation has no cost in energy 
in the fluid membranes. 

It was reported that a collapsing transition and a tran- 
sition of surface fluctuations occur in the conventional 
fixed connectivity surface models, where no in-plane shear 
energy is included in the Hamiltonians [TU^. These two 
transitions are identical to the so-called crumpling transi- 
tion, which has been extensively studied so far pn[T^[|T51 

mm- 

However, it is unclear whether the in-plane shear elas- 
ticity is negligible to the bending elasticity in membranes 
in the gel phase or in membranes supported by cytoskele- 
tons. In fact, red blood cells are known to have non- 
negligible in-plane shear elasticity [TBIITT]. 



Therefore, it is natural to ask whether the in-plane 
shear elasticity influences the transitions, which can be 
seen in the models defined without the in-plane shear en- 
ergy. In this article, we study two types of meshwork mod- 
els in |18j and [19j by including an in-plane shear energy 
in the Hamiltonian in order to see possible influences of 
the in-plane elasticity on the phase transitions. 

The first meshwork model, studied in [T^, is distin- 
guished from the conventional surface models because of 
the difference in the phase structures of them. The mesh- 
work model undergoes a first-order collapsing transition 
between the collapsed phase and the smooth phase and a 
continuous transition of surface fluctuations at the same 
transition point, while both of the transitions arc of first- 
order in the conventional models. The transition of surface 
fluctuations distinguishes the meshwork model in |18j from 
the conventional surface models such as those in [10] . 

The second meshwork model, studied in [TH], is ob- 
tained from the model in [18] by replacing the Gaussian 
bond potential with the Nambu-Goto potential, which is 
given by the summation over the area of triangles. The 
Nambu-Goto surface model is well-known as an ill-defined 
model in the sense that the model has no smooth surface in 
the whole range of the bending energy b when the Hamil- 
tonian includes the conventional two-dimensional bend- 
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ing energy of the type 1 — n • n as the curvature energy, 
where n is a unit normal vector of a triangle [2^. The 
well-definedness of the model in [19] is due to the one- 
dimensional bending energy, and a variety of shapes makes 
the model very different from the conventional ones in 
[10] , which have only the smooth phase and the collapsed 
phase. The models in [T81IT9] are also different form the 
surface models with cytoskeletal structures in .21.22J , be- 
cause the length L between the junctions of the models in 
[TSIITQ] is L = 1 in the unit of bond length while that of 
the models in [^TU^ is L> 1. 

We should note that the linear bending energy as- 
sumed in the meshwork models in [18j and |19j produces 
no resistance force against the in-plane deformations of the 
surface; the junctions of the meshwork play only a role for 
binding the one-dimensional skeletons. This makes us to 
expect that the in-plane energy, as an additional energy 
term in the Hamiltonian, has a non-trivial influence on 
the phase transitions of the models. For this reason, we 
study in this paper the meshwork models in [18] and [19j 
by including the in-plane energy in the Hamiltonians. 

Our main conclusion in this paper is that the in-plane 
shear energy enhances the surface fluctuations. We should 
comment on why the in-plane energy makes the transi- 
tions strong. The reason seems that the energy of thermal 
fluctuations is accumulated mainly in the bending defor- 
mation energy of the surface if the surface has large in- 
plane shear resistance. Since the surface deformation can 
be divided into the bending deformation and the in-plane 
deformation, then the bending deformation is enhanced if 
the in-plane deformation is suppressed. 

2 Models and Monte Carlo technique 

The triangulated lattices, on which two types of models 
are defined, are identical to those in [l8llT9] . By splitting 
the icosahedron, we have a triangulated spherical mesh- 
work of size N = 10-^^-1-2, which is the total number of 
vertices. The symbol £ in is the number of partitions 
of an edge of the icosahedron. The co-ordination number 
g, which is the total number of bonds emanating from a 
vertex, is g = 5 on 12 vertices, which correspond to those 
of the icosahedron, and g = 6 on the remaining A^ — 12 
vertices. 

The Hamiltonian S is given by a linear combination 
of the bond potential 5i, the one-dimensional bending 
energy 5*2, and the in-plane shear energy S^, which are 
defined by 

c _/'5'i,G (model 1), a ^^ + +\ 

'^'"l^i.NG (model 2), '^'2-2.(l-t.-t,), 

53 - 5] [1 - cos(5, - n/3)] . (1) 

i 

Si^G and 5*1, NG denote the Gaussian bond potential and 
the Nambu-Goto potential, which are respectively given 

by 

Si,G^Y.(^^-^o)\ ^l.NG^^A^. (2) 



The symbol Xi in 5i.g denotes the three-dimensional po- 
sition of the vertex i, and in 'S'l.G denotes the sum 
over bonds (ij), which connect the vertices i and j. 
in 5*1, NG denotes the area of the triangle A. We call the 
model defined by 5*1 = Si^q, 5*2, and ^3 as model 1, and 
the one defined by Si =5'i,ng, 'S'2, and 53 as model 2. 




(a) (b) (c) 

Fig. 1. (a) Unit tangential vectors ti and t2 at a vertex of 
coordination number q = 6, and (b) those at a vertex of coordi- 
nation number q = 5, and (c) three internal angles 5i{i = 1, 2, 3) 
in the definition of 5*3. The vectors in (a) and those in (b) de- 
fine the bending energy 1 — ti ■ t2 and [1— ti ■ (t2 +t3)/2]/2, 
respectively. 

In S'2 of Eq. ([T]), ti is a unit tangential vector of the 
bond i. The definition of in S'2 is identical to that in 

[TSlfTO] and is summarized as follows: 1— ti •t2 is included 
in S2 with the weight of 1 at the vertex of g = 6 shown in 
Fig. (TJa), while 1— ti • (t2-l-t3)/2 is included in S'2 with the 
weight of 1/2 at the vertex of q = 5 shown in Fig. (IJb) . As 
a result, we have ^ = ^b, where Nb = 3N — 6 is the 

total number of bonds. 

5i in S'3 of Eq. ([T]) denotes one of the three internal 
angles of a triangle such as shown in Fig.[ljc). Thus, the 
summation ^ - in S3 gives ^ - 1 = 3Nt, where A'^r(= 2A'^-4) 
is the total number of triangles. Although the angles 6i{i = 
1,2,3) are linearly dependent to each other, we use three 
of them in S3 because of its non-linear definition with 
respect to the angle 6. 

The models; model 1 and model 2, are defined by the 
following partition function: 

Z= l[dX,cxp[-S{X)], (3) 

1=1 

S{X)^ Si+bS2 + aS3, 

Si — Si^G (model 1), or Si ~ Si^ng (model 2), 

where b is the bending rigidity, and a is the parameter 
denoting the in-plane rigidity of the surface. In the limit 
of a ^ cxD, we expect that the surface is composed of only 
regular triangles. On the contrary, the meshwork models in 
[l8] and [l9j are restored in the limit of a — > 0. The prime 
in j' denotes that the integration is performed under the 
condition that the center of mass of the surface is fixed. 

We note that both b and a in Eq. ([3]) are the micro- 
scopic quantities and not always identical to the macro- 
scopic ones. So it is unclear to what value a should be fixed 
in the simulations to see influence of the in-pl|ine elasticity 
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on the transitions. However, as mentioned in the introduc- 
tion, we know that there is a membrane whose in-plane 
shear modulus is not negligible. We should remind our- 
selves of that the macroscopic bending modulus and the 
macroscopic shear modulus are of the same order in the 
red cells p^ll7| . Therefore, we assume the value of a in 
the simulations as a = 1, which is of the same order as 
that of b = bc the collapsing transition point. 

The surface has no in-plane shear resistance against 
the in-plane shear deformations in the case of a —s- in 
both models. The variation of the internal angles S in Fig. 
[T] changes neither Si nor 5*2. It is easy to understand that 
only S3 can reflect the variation of 6. In fact, thin and 
oblong triangles form the surfaces in the planar phase and 
in the linear phase in the model of [in]. The vertices of 
the surface in the models of [TSlfTO] have only a role for 
binding the one-dimensional skeletons. To the contrary, 
the surface has the in-plane shear resistance in the case 
that the junctions are elastic or rigid [22] . 



The dynamical variables X in Z are integrated over 
by the canonical Metropolis Monte Carlo technique. A 
random shift X X' = X+SX is accepted with the prob- 
ability Min[l,exp(-55)], where 5S' = 5'(new)-S'(old). The 
vector SX is chosen randomly in a sphere, whose radius is 
fixed at the beginning of the simulations for maintaining 
about 50% acceptance rate. 



Total number of Monte Carlo sweeps (MCS) after the 
thermalization MCS at the transition region of model 1 
is about (2 - 2.5) x 10^ for the N = 2562 surface, (4 ~ 
5) X 10^ for the N = 4842 surface, (6 - 7) x 10^ for the 
TV = 10242 surface, and (8 - 9) x 10** for the N = 16812 
surface. Relatively smaller number of MCS is performed 
at non-transition regions of b. The total number of MCS 
for model 2 is almost identical to or slightly smaller than 
that of model 1. 



The standard error Sq of the quantity Q is defined 
by the so-called the binning analysis: Let A'tot be the to- 
tal number of data in the sequence {Q}i{i = 1, A'tot), 
where each Qi is obtained at every 1000 MCS. The total 
number of MCS performed is thus given by lOOOA'^tot- The 
series {Q}i is split into rib sub-series {Q}i{I = l,...,ni,), 
and Qi denotes the mean value of the I-th sub series. 
The total number of MCS ntot in a sub-series is given by 
"■tot = lOOOiVtot/'^fc. Let Q denotes the mean value of {Q}i; 
Q is also the mean value of {Q}i- Then, 5q of {Q}i is de- 
fined by 6q = ^y^'i^i{Qi"Qy/nb, and Sq is expected to 
decrease with increasing ni, if ntot is sufficiently large such 
that {Q}i are statistically independent. The total number 
of MCS ntot in a series after the thermalization MCS is 
assumed as follows: ntot = 1 x 10'' for the N — 2562 sur- 
face and ntot = 2 x 10^ for the N = 4842, N = 10242, and 
A'' =16812 surfaces. The assumed value ntot for the large 
surfaces is insufficient and then, Sq remains large. 



3 Results 

3.1 Gaussian bond potential model 

In this subsection, the numerical data of model 1 will be 
presented. First we show in Fig.[2l^a) the mean square size 
X^ defined by 



X' 



-E 



X, - X] 



X 



-E 



X, 



(4) 



where X is the center of mass of the surface. The solid 
lines are drawn by the multihistogram reweighting tech- 
nique [23 ■ Rapid change of X"^ against b can be seen in 
the figure as iV increases, and this indicates the existence 
of the collapsing transition between the smooth phase and 
the collapsed phase. The transition appears to be discon- 
tinuous, because X"^ seems to change discontinuously on 
the largest surface. 




0.565 0.57 , 0.575 0.565 0.57 , 0.575 

(a) ^ (b) " (C) 



5000 10000 
N 



Fig. 2. (a) The mean square size versus h of model 1, (b) 
the variance C^a of X'^ versus 6, and (c) log-log plots of the 
peak C^^i^ against N . The solid lines in (a) and (b) are drawn 
by the multihistogram reweighting technique i23.. The straight 
line in (c) is obtained by fitting the data to Eq. ((6} . The critical 
exponent a = 1.36(1), which is given by the slope of the fitted 
line, indicates that the collapsing transition is of first-order. 



The variance of X^ is defined by 

cx^ = ^{{x^-m)"). 



(5) 



and this reflects how large the fluctuation of X"^ is. We 
show versus 6 in Fig. [IJb) in order to see the order 
of the transition more clearly. The solid lines are drawn 
also by the multihistogram reweighting technique. Sharp 
peaks seen on imply a phase transition. Therefore, 
we plot the peak values 6*^2'^, which are obtained by the 
multihistogram reweighting technique, in a log-log scale 
against X in Fig. [He). The error bars in Figs. [Ha) de- 
note the standard errors (5x2, while the errors in Fig.[Hb) 
are given by 5c^2 / \f^\ ^c^i is too large to show, hence 
we divide it by y/n\,. The error bars in Fig. [^c) denote 
the standard errors ^c^2 obtained by the multihistohgram 
technique. The fitting of data are performed by using 
those errors in Fig. [Ifc). 

The straight line in Fig. [DJc) is drawn by the power 
law fitting of the data, and we have a critical exponent cr 
such that 



C^f cx X" 



cr = 1.36 ± 0.01. 



(6) 
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Thus, the result confirms that the transition is of first- 
order from the finite-size scafing (FSS) theory [241125]. The 
exponent a is larger than one, i.e., cr > 1, however, the 
anomalous property C^f^ oo{N ^ cx)) is evidently 
seen. 
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Fig. 3. The variation of versus MCS at the transition 
region of (a), (b), (c) the A'' = 4842 surface, and (d), (e), (f) 
the iV = 10242 surface, and (g), (h), (i) the iV = 16812 surface. 
The dashed lines drawn in the figures denote the lower and the 
upper bounds for the evaluation of the mean values of in 
the collapsed phase and in the smooth phase. 

Figures [3la)-[3l^i) show the variation of against 
MCS of model 1 at the transition region of the surfaces 
A^ = 4842, Af= 10242, and iV= 16812. We see in Fig.He) a 
jump of X^ from the collapsed phase to the smooth phase 
and vise versa. If we compute the mean values of X^ in 
both phases separately, the Hausdorff dimension H can be 
evaluated by 



X' 



(N oo). 



(7) 



^2 

'min 

ma^x and Xj^jjj are shown in 



The dashed lines drawn horizontally in the figures denote 
the upper and the lower bounds X^^^, X^^^^ for evaluating 
the mean values of X^ in each phase such that X^;„ < 
X^ < X^^^. The values of 
Table [H 

The mean values of X^ computed by using the upper 
and the lower bounds in Table [T] are shown in Fig. Ufa) 
in a log-log scale. The error bars in Fig. HJa) denote the 
standard deviation obtained from the mean values of X^ 
and the series {X^}i with the condition imposed by X^^^ 
and X^-^. The straight lines drawn in Fig. IDJa) are ob- 
tained by fitting the data to the scaling relation of Eq. 
([7]). In the collapsed phase, the largest three X'^ are used 
in the fitting. Thus, we have 

(smooth), (8) 
(model 1). 

In the smooth phase, we have an expected result, which 
is close to the topological dimension H = 2, while we have 



= 2.22 ±0.15, 
4.1 ±2.8, (collapsed), 



Table 1. The lower bound and the upper bound X, 

for the mean value X^(col) in the collapsed state of model 1, 
and the lower bound X^^^f^° and the upper bound for 
the mean value X^{smo) in the smooth state. 



N 


b 


Y-2 col 
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0.57 
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0.566 


17 
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0.568 
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124 


16812 


0.5646 


22 


100 






16812 


0.566 






120 


195 



an unphysical H{> 3) in the collapsed phase. In the case 
of a = in [T^ , we have a physical H even in the collapsed 
phase. Therefore, the result H'^°^{> 3) in Eq. ([8]) allows us 
to understand that the in-pane share elasticity strengthens 
the transition at a = 0. Our model is not self-avoiding, 
consequently, the size of the surface in the collapsed phase 
shrinks if the transition is strong. As a consequence, the 
Hausdorff dimension becomes unphysical in the collapsed 
phase. 
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a=l 
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V:N=2562 
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Fig. 4. (a) Log-log plots of X^ vs. A'' obtained in the smooth 
phase and in the collapsed phase at the transition point of 
model 1, (b) the Gaussian bond potential Si^a/N vs. b. 



FigureHIb) shows the Gaussian bond potential Si^g/N 
vs. b. From the scale invariant property of the partition 
function, the relation Si,g/N — 3/2 should hold in the 
whole region of b and a. In our simulations, the expected 
relation is satisfied as we see in Fig. lUJb). Therefore, we 
consider that the simulations were correctly performed. 

Here we should note on the reason why a in Eq. ([6]) is 
larger than <j=l. The reason seems to be the low statis- 
tics of the simulations on the large sized surfaces. We see 
from Fig. [3l^e) that ntot should be ntot > 4 x 10^ or more 
even on the iV= 10242 surface. This implies that the high 
statistics simulations for surface models are not so easy 
even on the surface of size N = 10242. Nevertheless, we 
must emphasize that the first-order nature of the transi- 
tion is conclusive from the numerical results in this paper. 

The bending energy S2/NB, which is defined in Eq. 
(JlJ, is plotted in Fig. [51(a), where Nb is the total number 
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Fig. 5. (a) The one-dimensional bending energy S2/NB versus 
b of model 1, (b) the specific heat Cs2 versus b, and (c) log- 
log plots of the peak C^^^ against A''. The straight line in (c) 
is obtained by fitting the largest three data to Eq. (|10p . The 
critical exponent f = 1.13(8) indicates that the transition of 
surface fluctuations is of first-order. 

of bonds. We also see a rapid change in S2/NB just like in 
in Fig. [^^a). This suggests a discontinuous transition 

of surface fluctuations. 

The specific heat Cs^ for S2 is defined by 



S2 



N 



(9) 



and is plotted in Fig. [5jb) against b. An anomalous peak 
Cg^^^ can also be expected in Csa just as in Cx^ , and 
the peaks are apparently seen in Fig. [Sl^b). Figure (Sjc) 
shows a log-log plot of C^;^^ against N. The straight line 
is obtained by fitting the largest three data, and we have 
a critical exponent such that 

= 1.13 ±0.08, (model 1). (10) 



The result indicates that the transition of surface fluctua- 
tions is of first-order, because v almost satisfies i' = l. This 
is a remarkable result distinguishing model 1 in this paper 
from the model in [18] , where the transition of surface fluc- 
tuations is reported to be of second-order. Therefore, we 
understand that the in-plane shear elasticity strengthens 
the transition of surface fluctuations in the surface model 
with one-dimensional bending energy. 




Fig. 6. (a) The two-dimensional bending energy Sa/Nb versus 
b of model 1, (b) the variance Cs^ of S4 versus b, and (c) log-log 
plots of the peak C^^^ against A''. The straight line in (c) is 
obtained by fltting the largest three data to Cg^^ ~ A'^'^, and 
the result v = 1.10(8) indicates that the transition of surface 
fluctuations is of flrst-order. 



In order to confirm more convincingly this fact that the 
model undergoes a first-order transition of surface fluctu- 



ations, we plot in Fig. [6l[a) the two-dimensional bending 
energy defined by 



SA = ^{l-n^- rij) , 



(11) 



where is a unit normal vector of the triangle i. Although 
the bending energy 6*4 is not included in the Hamiltonian, 
it represents how large the surface fluctuates. The variance 
Cs, = {l/N){ {S4,-{Si)f) of 5*4 is plotted in Fig. Eljb), 
and we see that Cs^ has an expected anomalous peak, 
which reflects a discontinuous nature of S'4. Figure [6l|c) 
is a log-log plot of the peak value C^^^ against N, where 
the straight line is drawn by the fitting the data according 
to the scaling relation C™^^ ^ N'^, and we have a critical 
exponent fj, = 1.10 ± 0.08. Then, we again find that the 
transition of surface fluctuations is of first-order from the 
obtained exponent and the FSS theory. 



S3/(3Nt) 




0.565 0.57 , 0.575 
(a) (b) 



Fig. 7. (a) The in-plane shear energy S3/{3Nt) versus b, (b) 
the variance Csg of S3 versus b, and (c) log-log plots of the 
peak C^^^ against A''. The straight line in (c) is obtained by 
fitting the largest three data to C^^"" ~ A^*. 



Finally in this subsection, we show the in-plane energy 
S3/{3Nt) versus b, the variance Css of S3 versus 5, and 
the log-log plot of the peak value Cg^^^ against N, in Fig. 
[7l[a), Fig.[7Ub), and Fig.[7fc), respectively. An anomalous 
peak can also be seen in Cs^ ■ The scaling behavior can be 
seen in Fig.[7]Jc), and we have 



C 



S3 



(5 = 0.76 ±0.41. 



(12) 



This implies that the in-plane order- disorder transition is 
of first-order because 5 can be seen as S — 1 within the 
error, however this conclusion is less accurate because of 
the large errors. 

The in-plane order-disorder transition is expected to 
be discontinuous at least when a = 1 in Eq. ^ , and both 
of the collapsing transition and the transition of surface 
fluctuations are of first-order. Moreover, we expect that 
the order of the in-plane order-disorder transition changes 
according to the value a. The transition is expected to 
change from the discontinuous one to a continuous one at 
certain value of a at least in < a < 1. 



3.2 Nambu-Goto potential model 

The numerical data obtained from model 2 are presented 
in this subsection. The presentation is almost parallel to 
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the one of model 1 in the previous subsection. The pa- 
rameter a in Eq. Q is fixed to a = 1 just as in model 
1. 



"(a) |m=4S42. b=Q.64.c^l| 
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Fig. 8. (a) The mean square size versus b of model 2, (b) 
the variance (7^2 of versus b, and (c) log-log plots of the 
peak against A*'. The solid lines in (a) and (b) are drawn 

by the multihistogram reweighting technique [23]. The straight 
line in (c) is obtained by fitting the data to Eq. The critical 
exponent a = 1.11(9) indicates that the collapsing transition is 
of first-order. 



Figure [5Ja) shows the mean square size defined in 
Eq. ([4]). The solid curves are also drawn by the multihis- 
togram technique like those in the figures in the previous 
section. The variance Cx2 is shown in Fig. [Sfb), and the 
peak values Cx2^ are plotted in Fig. [8jc) against in a 
log-log scale. The straight line is drawn by fitting all of 
the four data to Eq. We have the critical exponent 
(T= 1.11±0.09, which implies a first-order collapsing tran- 
sition. 

We should remark that the in-plane energy 53 de- 
creases the multitude of phases of the model in [12] , and 
only two phases; the smooth and the collapsed phases, sur- 
vive. These two phases always seen in the standard surface 
models such as those of 10 and in the mesh work model 
of [18]. The reason for the shrinkage in the multitude of 
phases in the model of [T^] is intuitively understood; in 
fact, the linear and the planar surfaces are constructed 
not only from regular triangles but also from oblong tri- 
angles. To the contrary, oblong triangles are prohibited to 
occur in model 2 of this paper due to the presence of ^3. 

Figures [ni[a)-(5Ki) show the variation of against 
MCS of model 2, where N = 4842, N = 10242, and N = 
16812. It is clear that two different states coexist; one is 
characterized by the large X^ and the other by the small 
X'^, which respectively correspond to the smooth phase 
and the collapsed state. The transition, which separates 
one state from the other, can be called the collapsing tran- 
sition just as in model 1. 

In order to see the first-order nature of the collapsing 
transition more clearly, we show in Figs. fToral-fTOTc) the 
normalized distribution (or histogram) h{X^) of X^ at 
the transition point of the surfaces of = 4842, N = 
10242, and iV = 16812. The histograms in Figs.[TOi:a)-[TO|;c) 
correspond to the variations of X^ in Figs.dj^e), anddjjh), 
respectively. The double peak structure can be seen in 
h{X'^) of the AT = 10242 surface and the iV= 16812 surface. 
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Fig. 9. The variation of versus MCS of model 2 at the 
transition region of (a), (b), (c) the N = 4842 surface, and 
(d), (e), (f) the TV = 10242 surface, and (g), (h), (i) the TV = 
16812 surface. The dashed lines drawn on the figures denote 
the lower and the upper bounds for the mean values of X^ in 
the collapsed phase and in the smooth phase. 
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(a) 

Fig. 10. The distribution (or histogram) 
transition point of the surface of size (a' 
10242, and (c) Ar = 16812 of model 2. 



h{X'-^) of X^ at the 
TV = 4842, (b) TV = 



and this clearly indicates that the collapsing transition is 
of first-order. 

It is interesting to see whether the collapsing transi- 
tion is physical or not, i.e., whether or not the Hausdorff 
dimension H is less than the physical bound in the col- 
lapsed phase. The model is allowed to self-intersect like 
model 1 in the previous section, however, the transition 
is considered to be physical if both of H in the smooth 
phase and the collapsed phase are less than the physical 
bound; iJ < 3. In Table[2l we show the lower and the upper 
bounds for the mean value X^, from which we have H by 
using Eq. ([7]). The values shown in Tableware indicated 
by the horizontal dashed lines in Figs. [9l[a)-[9](i) . 

Figure fTlT a) shows log-log plots of the mean values 
X^ against TV, where X^ were obtained from the variation 
of X^ shown in Figs. [9Ka)-[9]^i) by using X^^^^ and X^-^ 
in Table [2j The straight lines were drawn by fitting the 
largest three data to Eq. ([7]) both in the smooth phase 
and in the collapsed phase. We have 



2.13 ±0.27, (smooth). 



(13) 
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Table 2. The lower bound Xj^™' and the upper bound X^S' 
for the mean value X^(col) in the collapsed state of model 2, 
and the lower bound X^f™° and the upper bound X^^° for 
the mean value X^(smo) in the smooth state. These values 
correspond to the horizontal dashed lines in Figs. [9l[a)-(9]^i). 



N 


b 
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■y2 smo 
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4842 


0.64 


50 


160 






4842 


0.66 
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10242 


0.66 






320 


530 


16812 


0.65 


200 


480 






16812 


0.654 






480 


800 




5000 10000 - ^ 0.65 0.67 0.69 

(a) N (b) ^ 

Fig. 11. (a) Log-log plots of vs. A'^ obtained in the smooth 
phase and in the collapsed phase at the transition point, (b) 
the Nambu-Goto potential Si^-^a/N vs. b. The fittings in (a) 
were done by using the largest three data both in the collapsed 
phase and in the smooth phase. 

H''°^ = 1.95 ± 0.44, (collapsed), (model 2). 

Thus, the collapsing transition of model 2 is considered to 
be physical, because 7?'^°' is less than the physical bound. 
This point distinguishes model 2 from model 1, whose col- 
lapsing transition is considered as unphysical at least when 
a = l as demonstrated in the previous subsection. 

The Nambu-Goto potential ^i^ng/^ is plotted in Fig. 
[TlTb). We find that the expected relation S'i,ng/^ = 3/2 
is satisfied. This implies that model 2 is well-defined and 
moreover that the simulations of model 2 are successfully 
performed. 

Finally in this subsection, we plot S2/NB vs. b in Fig. 
[T^a) and the specific heat vs. b in Fig. [TWb). The 
peaks Cg]^^ of the specific heat seen in Fig. [T2r b) imply 
that model 2 undergoes a phase transition of surface fluc- 
tuations. The log-log plot of Cg]^^ vs. N is shown in Fig. 
[T2lc). The exponent v defined by Eq. pO|) is obtained as 
follows: 

= 0.45 ± 0.06, (model 2). (14) 

This value of v implies that the transition of surface fluc- 
tuations is of second-order at a^l. However, a possibility 
of the first-order transition is not completely eliminated. 
Large scale simulations are still necessary to clarify this 
point. 




0.65 0.66 . 0.65 0.66 . 500010000 

(a) (b) (C) N 



Fig. 12. (a) The one- dimensional bending energy S2/NB vs. 
b of model 2, (b) the specific heat Cs^ versus b, and (c) log-log 
plots of the peak C™" against N . The straight line in (c) is 
obtained by fitting the data to Eq. (|10[) . The critical exponent 
= 0.45(6) indicates that the transition of surface fluctuations 
is of second-order. 

We performed the simulations for model 2 with a = 0.5 
by using the same sized surfaces as those for the case a = 1 
presented above. The phase structure is almost identical 
to the case a = l except the strength of the transitions. 
The collapsing transition is considered to be weakened but 
it still remains discontinuous at a = 0.5, because the finite 
size-scaling analysis of the peak values of C^f^ indicates 
a first-order collapsing transition. However, the transition 
of surface fluctuations changes to a higher-order one and 
almost disappears. Thus, we find that the strength of the 
transitions changes depending on the value of a in model 
2 as well as in model 1. 



4 Summary and conclusions 

In this article, we have investigated how the in-plane shear 
elasticity influences the transitions observed in two-types 
of the meshwork models in [TB] and pjl. The models are 
described as follows: The surface shape of the first model 
in [18] is maintained by a one-dimensional bending en- 
ergy and by the Gaussian bond potential, while the sur- 
face shape of the second model in [l9j is maintained by 
the one-dimensional bending energy and by the Nambu- 
Goto potential. One of the transitions in the model of [TH] 
is called the collapsing transition, which is of first-order 
and separates the smooth phase from the collapsed phase 
at finite bending rigidity. The other is a continuous one 
called the transition of surface fluctuations. We should re- 
mark that these two transitions are of discontinuous in the 
conventional surface model with two-dimensional bending 
energy jJO] . The second model studied in [19^ has a variety 
of phases including the smooth, the planar, and the linear 
phases in contrast to the first model in [18]. 

The reason why we study in this paper the meshwork 
models of [IB] and [T^] is that the models have no in- 
plane resistance force against the in-plane deformations. 
The models have a resistance only against the bending 
deformation and the tensile deformation. The linear sur- 
face of the model in [19] consists of oblong surfaces, and 
the planar surface consists of both regular triangles and 
oblong ones. No energy cost is necessary for the in-plane 
deformations in both of the models in [TBlfTO] . For this 
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reason, we expect that the in-plane energy makes a non- 
trivial effect on the transitions. We call the first and the 
second models with the in-plane energy as model 1 and 
model 2, respectively. The value of a the coefficient of the 
in-plane energy was assumed as a = 1 in the simulations 
in both models. 

Our numerical results of model 1 show that the tran- 
sition of surface fluctuations is of first-order, which was 
reported to be of second-order in the model without the 
in-plane shear energy just stated as above. The collaps- 
ing transition is also strengthened, although the order of 
the transition is identical to the case without the in-plane 
energy. Therefore, we conclude that the in-plane shear en- 
ergy can strengthen the transitions observed in the mesh- 
work model with the Gaussian bond potential. 

Moreover, we find in model 1 that the in-plane order- 
disorder transition is of first-order; the in-plane energy 
discontinuously changes against the bending rigidity at 
the transition point, where both of the one-dimensional 
and the two-dimensional bending energies also discontinu- 
ously change. The first-order nature of the in-plane order- 
disorder transition is not so accurate because of the large 
errors in the critical exponent for the variance of the in- 
plane energy. 

In the case of model 2, the variety of phases seen in 
the model in [19j disappears, and the phase structure is 
almost identical to that of model 1. The reason of this is 
because oblong triangles are suppressed due to the pres- 
ence of the in-plane energy. The collapsing transition be- 
tween the smooth phase and the collapsed phase is of first- 
order. Moreover, this transition is considered as physical, 
because the Hausdorff dimension H is less than the phys- 
ical bound; H < 3, even in the collapsed phase close to 
the transition point. The transition of surface fluctuations 
of model 2 is considered to be continuous. The finite-size 
scaling analysis of the specific heat of the one-dimensional 
bending energy supports this conclusion. 

We should note that the strength of the transitions 
changes depending on the value of a. One can also ex- 
pect that the unphysical collapsing transition of model 1 
changes to the physical one at some values of a in the 
range < a < 1, because the strength of the transitions 
weakens with decreasing a. We know that the collapsing 
transition of model 1 is physical in the limit of a — > 
[TB]. It is also possible that the transition of surface fluc- 
tuations in model 2 turns to be a discontinuous one with 
increasing a. It is interesting to study the conventional 
curvature surface model in |10| by including the in-plane 
shear energy in the Hamiltonian. 
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